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Abstract. Infrared spectra of a number of evolved O-rich 
stars have been obtained with the Short- and Long- Wave- 
length spectrometers on board the Infrared Space Obser- 
vatory. The very broad wavelength coverage (2.4-197/im) 
obtained by combining observations made with the two 
spectrometers includes practically all of the flux emitted 
by the sources, and allows us to determine the emission 
and absorption features of the dense circumstellar dust 
shells. Agreement between the fluxes obtained by the two 
instruments is generally very good; the largest discrepan- 
cies are probably due to source variability. Our sample of 
oxygen-rich AGB stars exhibits a wealth of spectral fea- 
tures due to crystalline silicates and crystalline water ice 
in emission and absorption. In this study a qualitative 
overview of all features due to crystalline silicates and wa- 
ter ice in these high mass loss rate objects is presented. It 
seems that there is a certain onset value for the mass loss 
rate above which these features appear in the spectrum. 
Moreover, crystalline silicate emission features have been 
detected for the first time at wavelengths where the amor- 
phous silicates are still in absorption, implying different 
spatial distributions for the two materials. A spherically 
symmetric and an axi-symmetric geometry are proposed. 

Key words: Stars: AGB and post- AGB - circumstellar 
matter - Stars: mass-loss - Infrared: stars 
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1. Introduction 

The post-main sequence evolution of stars of low or inter- 
mediate mass takes these stars on to the Asymptotic Giant 
Branch (AGB), where they lose mass at rates of 10~^- 
lO^''M0yr^^. In the circumstellar outflows molecules are 
formed and dust grains condense. The relative abundances 
of carbon and oxygen in the star determine the chemical 
composition of the gas and dust in the outflows. Oxygen- 
rich stars produce silicate dust and molecules such as 
H2O and OH. If the mass-loss rate is sufficiently high, the 
dust completely obscures the star at visible wavelengths, 
and the object is known as an OH/IR star because of 
its strong emission in the infrared (IR), produced by the 
dust grains, and in radio OH lines, due to maser action 
by OH molecules. See Habing (1996) for a detailed re- 
view of AGB and OH/IR stars. The optically-thick dust 
envelopes of OH/IR stars may be the result of a recent 
increase in mass-loss rate: the so-called 'superwind' phase 
(e.g. Justtanont et al. 1996b; Delfosse et al. 1997). Omont 
et al. (199C) detected the 43- and 60- /xm emission bands 
of water ice in the KAO spectra of a number of OH/IR 
stars, attributing the bands to the condensation of water 
molecules onto silicate grain cores in the dense outflows. 
An important result from the Infrared Space Obser- 
vatory (ISO) mission is the detection of emission from 
crystalline silicates in the far-IR spectra of many sources, 
including the circumstellar environments of young and 
evolved stars and solar-system comets (e.g. Waters et 
al. 1996, |1999 ). Crystalline silicate bands have been de- 
tected in OH/IR star spectra (Cami et al. 1998) but are 



not seen in O-rich AGB stars with low mass-loss rates, 
suggesting that the abundance of the crystalline materials 
is related to the density of the circumstell ar ma tter at the 
dust condensation radius (Waters et al. 1996). ISO has 
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also detected thermal emission and absorption by water, 
in both the gaseous and solid (ice) phases (e.g. Barlow 
199S), from 0-rich circumstellar environments. 



In this paper we present ISO spectra of seven well- 
known OH/IR stars covering a range of mass-loss rates. 
The spectrum of the archetypal Mir a variable, o Get, is 
presented for comparison. For most of our targets, the 
spectra cover the complete 2.4-197/Ltm spectral range of 
ISO. Sect. 2 of the paper describes the observations and 
data reduction. In Sect. 3, the spectra are presented and 
analyzed, with emphasis on the determination of the con- 
tinuum and the features due to ices and silicates. Con- 
cluding remarks are made in Sect. 4. 

2. Observations 

Seven of our eight sources were observed with both the 
SWS and LWS instruments, while of OH32. 8-0.3 only the 
SWS data were useful. Table |l| lists the sources observed, 
and the JD dates of the observations (henceforth we abbre- 
viate the OH/IR star designations to OH104.9, OH127.8 
etc). Some of our sources (e.g. OH104.9, OH26.5) were 
observed nearly contemporaneously with the two instru- 
ments, while for others, the two spectra were taken more 
than 100 days apart. In the past , modelling work (e.g. 
Lorenz- Martins & de Araujo 1997) has been hindered to 
some extent by the non-simultaneity of NIR photometry 
and 10-20 /im spectra (usually IRAS LRS data). The fact 
that the SWS spectrum covers both these wavelength re- 
gions at a single epoch will be useful for future modelling 
of these sources (Kemper et al, in preparation). 

2.1. SWS 

The 2.38-45.2 /im part of the spectrum was obtained us- 
ing the ISO Short Wavelength Spectrometer (SWS). A 
detailed description of the instrument can be found in de 



Graauw et al. (1996a). Our objects were observed in AOT 
1 mode, speed 2, except for Mira (speed 3) and OH104.9 
(speed 1). The spectrum scanned with SWS contains 12 
subspectra, that each consist of two scans, one in the di- 
rection of decreasing wavelength ('up' scan) and one in the 
increasing wavelength direction ('down' scan). There are 
small regions of overlap in wavelength between the sub- 
spectra. Each sub-spectrum is recorded by 12 independent 
detectors. 

The data reduction was performed using the ESA SWS 
Interactive Analysis package (lA'^), together with the cali- 
bration files available in January 1999, equivalent to pipe- 
line version 6.0. We started from the Standard Processed 
Data (SPD) to determine the final spectrum, according to 
the steps described in this session. 

The observations suffer from severe memory effects in 
the 4.08-12.0 /zm and 29.0-45.2 /im wavelength regions. 
It is possible to correct for the memory effects for the 
individual detectors, using a combined dark-current and 



memory-effect subtraction method. This method was ap- 
plied assuming that the flux levels in these wavelength 
regions are very high, and treating the memory effect as 
giving an additive contribution to the observed signal. We 
also assumed that the spurious signal from the memory 
effect reaches a certain saturation value very quickly af- 
ter the start of the up scan, and then remains constant 
throughout the rest of the up scan and the entire down 
scan. This memory saturation value is measured imme- 
diately after the down scan is ended, and is subtracted 
from the up and down scan measurements. The spectral 
shape of the memory-corrected down scan is now correct; 
the error in the spectral shape of the up scan is corrected 
by fitting a polynomial to the up scan and adjusting this 
fit to the down scan, without changing the detailed struc- 
ture of the spectrum. The order of the applied polynomial 
fit differs per subband, but is chosen to be in agreement 
with the spectral shape in that band. For band 1 and 4 we 
mostly used polynomials of order 1 or 2, for band 2a, 2b 
and 3 we predominantly used order 2 or 3, and for band 2c 
higher order polynomials (up to order 10) were required 
to adjust the up scan to the down scan. 

The spectra of some of our objects showed fringes in 
the 12.0-29.0 /itm wavelength region. This was corrected 
using the defringe procedures of lA"^ . 

Glitches caused by particle hits on the detector were 
removed by hand. Glitches can be easily recognized: they 
start with a sudden increase in flux level, followed by a 
tail which decreases exponentially with time. Any given 
glitch affects only one of the two scans. 

The data were further analyzed by shifting all spectra 
of the separate detectors to a mean value, followed by 
sigma clipping and rebinning to a resolution of A/AA = 
600, which is reasonable for AOT 1 speed 2 observations. 

2.2. LWS 

We obtained 43-197 /xm grating spectra using the LWS in- 
strument. Details of the instrum ent a nd its performance 
can be found in Glegg et al. ( 1996| ) and Swinyard et 
al. ( 1996 ) respectively. The resolution element was 0.3 /xm 
for the the short-wavelength detectors (A < 93 ^m) and 
0.6 ^m for the long- wavelength detectors (A > 80 /im). 
Four samples were taken per resolution element. Between 
6 and 26 fast grating scans were made of each target, de- 
pending on source brightness and scheduling constraints. 
Each scan consisted of a single 0.5-sec integration per sam- 
ple. 

The data were reduced using the LWS off-line process- 
ing software (version 7.0), and then averaging the scans af- 
ter sigma-clipping to remove the discrepant points caused 
by cosmic-ray hits. 

For all our LWS sources, apart from Mira and WX 
Psc, the Galactic background FIR emission was strong, 
and off-source spectra were taken to enable the back- 
ground to be subtracted from the on-source spectrum. 
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Table 1. Observing details for our sources. Included are the revolution number and on-target time for the SWS and 
LWS observations, and the scale factor applied to the LWS observations to bring the spectra into agreement (see text). 



Star 


IRAS name 


Date SWS Obs 


Date LWS Obs 


LWS Factor 


SWS Time 


LWS Time 






JD 2450000+ 


JD 2450000+ 




(sec) 


(sec) 


Mira 


02168-0312 


489.47 


633.83 


1.5 


3454 


1928 


CRL 2199 


18333+0553 


749.17 


746.42 


0.97 


1912 


2228 


WXPsc 


01037+1219 


433.43 


614.52 


1.4 


1912 


4704 


OH104.9+2.4 


21177+5936 


321.29 


321.28 


1.0 


1140 


1268 


OH127.8+0.0 


01304+6211 


825.20 


651.02 


1.0 


1912 


1330 


OH26.5+0.6 


18348-0526 


368.16 


368.13 


1.15 


1912 


1268 


AFGL 5379 


17411-3154 


879.85 


507.23 


0.8 


1912 


1614 


OH32.8-0.3 


18498-0017 


358.42 


- 


- 


1912 


- 



Galactic background flux levels were only significant for 
A > 100 /im. The Galactic background emission is ex- 
tended compared to the LWS beam, and therefore gives 
rise to strong fringing in both the on- and off-source spec- 
tra. The background-subtracted spectra do not show fring- 
ing, indicating that the OH/IR stars are point-like to the 
LWS, as expected. 

After averaging and background subtraction (if neces- 
sary), each observation consisted of ten subspectra (one 
per detector) , which were rescaled by small factors to give 
the consistent fluxes in regions of overlap, and merged to 
give a final spectrum. 

2.3. Joining the SWS and LWS spectra 

One of the main goals of this article is to study the overall 
ISO spectra of the selected objects. Therefore, it is nec- 
essary to join the SWS and LWS spectra in such a way 
that the flux levels and slopes of the spectra agree for both 
LWS and SWS. Differences in the flux levels of the LWS 
and SWS spectra are mostly due to flux calibration uncer- 
tainties. Although the spectral shape is very reliable, the 
absolute flux calibration uncertainty is 30% for the SWS 
at 45 ^lui (Schaeidt et al. |1996| ), and 10- 15% for the LWS 
at the same wavelength (Swinyard et al. 1998 ). Therefore, 
differences between the flux levels of LWS and SWS which 
are smaller than 33% are acceptable within the limits of 
the combined error bars. 

The SWS and LWS spectra were scaled according to 
their fluxes in the overlap region. Generally this resulted 
in a shift of less than 20% (see Table |^). In the case of Mira 
and WX Psc, a much larger shift was required, presumably 
due to the large time interval between the SWS and LWS 
observations of these variable stars. 



3. Results 

The combined spectra are presented in Fig. |l|, in \F\ 
units. 

The spectra are ordered by increasing optical depth 
in the observed 10-/im silicate absorption, Tg, and hence 
are in approximate order of increasing mass-loss rate, as- 



Table 2. Dust mass loss rates and SED properties of our 
sample stars. Mass loss rates are from Schutte & Tielens 
( |1989D and Justtanont & Tielens ( |1992D . Tbb is the tem- 
perature of a blackbody which approximates the observed 
SED. Tg was measured directly from the spectra, using 
continuum fits derived according to the method described 
in paragraph 3.1, while Tmod is the silicate optical depth 



derived by Justtanont & Tielens (1992) from radiative 
transfer modelling. 



Star 



Mi (Mq yr-^) rBB(K) 



'^mod 



Mira 




1000 






CRL 2199 


1.8x10"'^ 


600 


1.2 


3.2 


WXPsc 


7.6x10"* 


550 


1.3 




OH104.9+2.4 


2.4x10-'^ 


400 


1.2 




OH127.8+0.0 


2.0x10-*^ 


370 


2.2 


16.0 


OH26.5+0.6 


1.2x10"'' 


350 


2.6 


19.6 


AFGL 5379 




300 


3.5 




OH32.8-0.3 


2.2x10-*^ 


280 


3.9 


19.0 



suming roughly similar luminosities. By modelling the in- 
frared excess emission of 0-rich AGB st ars, S chutte & 
Tielens (|l989| ) and Justtanont & Tielens ( |l992[ ) have de- 
termined the dust mass loss rates for several individual 
0-rich AGB stars. Their results are summarized in Ta- 
ble g, where our sources are listed in the same order as 
in Fig. |[ According to Table g, our sample is indeed or- 
dered with increasing mass loss rate, excluding WX Psc. 
The values of Tg measured from the spectra are the ap- 
parent optical depth compared to our continuum fit to 
the overall SED, or to an assumed silicate emission profile 
for CRL 2199 and WX Psc. They therefore represent only 
a part of the total lO-fiui silicate optical depth towards 
the sources. This is evident from Table ||, where the mea- 
sured Tg are quoted, along with the 10-/im optical depths 
derived for some of our sources by Justtanont & Tielens 



(1992), using radiative transfer modelling. A description 
of the determination of our continuum fit can be found in 
Sect. 3.1. 

All of our sources are dominated by continuum emis- 
sion from cool dust. The most striking feature in our spec- 
tra is the lO-fim silicate band, which appears in emis- 
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Fig. 1. Combined SWS and LWS spectral energy distributions of our programme stars. Flux units are log(Ai^A) in 
Wm^^. The spectra have been shifted along the ordinate to order them in terms of 10- /xm optical depth, which 
correlates roughly with mass-loss rate. The amount by which the spectra were shifted is given in parentheses. The 
continuum fit to OH104.9 is also included in the figure (dashed line). 



sion for Mira, partially self-absorbed for WX Psc and 
CRL 2199, and in strong absorption for the other sources. 
The overall shape of the observed spectral energy distribu- 
tion (SED) varies with Ts, the sources with high Tg showing 
'redder' SEDs. The SEDs can be roughly approximated us- 
ing blackbodies with temperatures ranging from 300 K for 



the most optically-thick sources, to 600 K for the partially 
self-absorbed sources (see Table ||). 

The 20-fim silicate feature also passes from emission 
to absorption going down the sequence in Fig. nl but none 
of our sources show it as a self-absorbed emission feature. 
Variations in shape as well as optical depth are evident 
in the 10- and 20-/im absorption features. Weaker fea- 
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tures beyond 20 /im can be discerned: these are features 
of water ice and crystalline silicates, and will be discussed 
below. Evident at short wavelengths are molecular absorp- 



tion bands (see Justtanont et al. 1996a for identifications 
for the supergiant source NML Cyg, which is not included 
here) and the 3.1-/zm H2O ice absorption feature. 

3.1. Determination of the continuum 

In order to determine the shape and relative strength 
of the emission and absorption features we will define a 
pseudo- continuum, which is assumed to represent fea- 
tureless thermal emission from the dust. The continuum- 
divided spectrum will provide us with information on the 
optical depth of the different species that are located out- 
side the continuum-producing region and the wavelength 
at which the material becomes optically thin. However, 
one has to be very careful with the physical meaning of the 
pseudo- continuum, since the observed continuum results 
from several wavelength-dependent dust emissivities, with 
a large range of dust temperatures, modified by strong 
optical depth effects, rather than being simply a superpo- 
sition of blackbody energy distributions corresponding to 
the physical temperature gradient. 

For the determination of the continuum we plotted the 
energy distribution log F\ (W m~^ /xm~^) against A (/xm). 
Plotting the data this way one can easily recognize the 
general shape of the continuum; the parts of the spectrum 
where the continuum is well defined, i.e. the long and short 
wavelength edges, are emphasized. This makes it easier to 
constrain the continuum in the wavelength regions where 
the strong spectral features are present. At A > 50 /xm, the 
dust is optically thin, and only weak emission features are 
present. At A < 7 /im, the radiation from the stellar pho- 
tosphere is partially (for the Miras) or completely (OH/IR 
stars) absorbed by molecular lines and the high dust opac- 
ity towards the central star. When a spline fit is performed 
in log F\ space, we can thus constrain the continuum by 
using the known long and short wavelength continuum 
contributions as reference points. Applying spline fitting 
with the same reference points in other parameter spaces 
(A log F\ space, F^, space, log F^, space etc.) shows that 
all the fits vary within 10% (in flux) with respect to the 
adopted pseudo-continuum in the 10-20 /im region. 

3. 2. Description of the spectrum 

Figs. 0-0 show the observed spectra after division by the 
adopted continuum. A wealth of interesting features can 
be recognized in the spectra. At the shortest wavelengths, 
from A = 2.38 - ^7 /im, see Fig. |[ the spectrum is dom- 
inated by molecular absorption bands and the effects of 
dust absorption. In Mira and the intermediate type stars 
(WX Psc and CRL 2199), strong absorption features due 
to a blend of H2O and CO are present from 2.38-'-^3.3 /im 
(2.38 — 3.8 /im in the case of WX Psc), see Fig. |[ Absorp- 



tion due to SiO is found in some objects around 4.0 /^m, 
most prominently in WX Psc and OH104.9. In most of our 
spectra we find clear absorption features of gaseous CO2 at 
4.3 /im. CO line absorption around 4.5 /im is found in all 
our spectra. Throughout the whole 2.38-7 /im wavelength 
region, H2O absorption lines are present; very strong lines 
of gaseous H2O are found at 6.60-6.63 /im. Yamamura et 



al. (1999) have performed a detailed study of the water 



features in the spectrum of Mira; Mira is also the only 
source in our sample that shows OH absorption lines in 
the 2.5-3.5 /xm region. In the high mass-loss rate objects, 
water ice features are present as well, at 3.1 and (possi- 
bly) 6.0 /im. The ice features will be discussed in the next 
section. 

The 7-25 /im region (Fig. 0) is dominated by the strong 
features of amorphous silicates. The 9.7-/im feature (rang- 
ing from 8 to 12-13 /zm) and the 18 /xm feature (ranging 
from 15 to 20-25 /im) are very strong bands that occur in 
emission or (self-)absorption. Crystalline silicate features 
also occur in the same wavelength ranges but are usually 
much weaker, so the actual shape of the observed silicate 
bands is due to a blend of crystalline and amorphous sili- 
cates. Absorption by molecular SiO can be seen at 8 /xm, 
on the wing of the silicate feature. Around 15 /im the 
spectrum of Mira exhibits some sharp emission lines due 
to CO2. 

Longwards of 25 /im (Fig. ^, most spectral features 
occur in emission. The only absorption feature is the OH 
pumping line at 34.6 /xm, which is detected in some of 
the sample stars. Groups of crystalline silicate emission 
features occur near 28, 33, and 43 /xm, the latter probably 
being a blend with the 43 /xm crystalline water ice feature, 
while a broad feature, which we ascribe to water ice, peaks 
around 62 /xm. Superposed on this is a sharp feature at 
69 /xm, due to forsterite. 

Several sources show emission features near 47.5 and 
52 /xm. The longitudinal optical band of water ice lies close 



to 52 /xm (Bertie & Whalley 1967), but in laboratory data 
it is only ever seen as a shoulder on the 43-/xm band (e.g. 
Smith et al. 1994), not a clearly-separated feature. The ob- 
served feature is therefore unlikely to be the longitudinal 



ice band. Malfait et al. (1999) suggested that montmoril 



lonite, which gave a good fit to the broad 100-/xm emission 
feature in the spectrum of HD 142527, is a possible carrier 
of the 47 and 50/xm features in this source. The planetary 
nebula NGC 6302 also shows the 47, 53 and 100 /xm bands 
(Lim et al., 1999). If the three features do have a common 
carrier, the apparent absence of the 100 /im band in the 
spectra presented here may be due to a difference in dust 
temperature. 

As well as broad emission bands, unresolved emission 
lines can be detected in the long-wavelength regions of 
the least noisy of our spectra, (e.g. WX Psc and AFGL 
5379) . Most of these lines are pure rotational lines of water 
vapour. The 157.7-/xm [C 11] line is also visible in most of 
our sources; this is the residual Galactic background [C 11] 
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Fig. 2. Continuum-divided spectra of our sources in the 2-7 fim region. The spectrum of OH32.8 is noisy at short 
wavelengths and has been omitted from this figure. The central wavelengths of absorption features due to H2O (gaseous 
and ice phases), and of gaseous CO, CO2 and SiO are indicated. Typically, the noise level (standard deviation) of the 
spectra is 2%, but is slightly larger for the short wavelengths of the high mass loss objects. 
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Fig. 3. Continuum-divided spectra of our sources in the 6-25/im region. The broad amorphous sihcate features at 
9.7 fim and 18 /im dominate this region of the spectra. Also visible are the SiO absorption at 8 fiia and the CO2 
emission lines around 15/xm, the latter only in the spectrum of Mira. The noise level (standard deviation) is typically 
1%, therefore we believe that most of the fine structure visible superposed on the broad silicate features is real. 
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Fig. 4. Continuum-divided spectra of our sources in the 20-198/im region. The central wavelengths of the broad 
emission features due to crystalline silicates and water ice are indicated. The noise level (standard deviation) of the 
data is about 2% up to 45 /im and about 4% at larger wavelengths. The spectrum of Mira shows poor agreement 
between the slopes of the SWS and LWS components, making continuum placement difficult in this wavelength region; 
it has therefore been omitted from this plot. There is no clear evidence for crystalline silicates or water ice in our Mira 
spectrum. 
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Table 3. Detections of water ice features in the spectra of 
our programme stars, ra.i and iVice are the optical depth 
of the 3.1-^m feature and the derived column density 



Source Absorption Emission 


TS.l 


iVieo 




3.1/im 


43pim 


62/im 




(lO'Vm-^) 


Mira 


N 


N 


N 


<0.02 


<3.4 


CRL 2199 


N 


N 


N 


<0.02 


<3.4 


WXPsc 


N 


Y 


N 


<0.02 


<3.4 


OH104.9 


N 


N 


? 


<0.05 


<8.4 


OH127.8 


Y 


Y 


Y 


0.33 


55 


OH26.5 


Y 


Y 


Y 


0.25 


42 


AFGL 5379 


Y 


Y 


Y 


0.67 


112 


OH32.8 


Yi 


Y 


- 


0.7^ 


120^ 


^ Roche & Aitken (1984 


) 








^ Meyer et al ( 


1998|) 





emission after subtraction of the off-source spectrum. In 
the case of AFGL 5379, the [C ii] line is seemingly in 
absorption: again, this is due to imperfect cancellation of 
the background emission. 

3.3. Water Ice 

Water ice is an important component of the solid-phase 
material in cool astronomical sources. Its spectrum shows 
bands at 3.1, 6.0, 11-12, 43 and 62 fim. The 3.1-fim 
stretching band is seen (always in absorption) in the spec- 
tra of many highly-embedded young stars (e.g. Whittet 
et al. 198S) and in some OH/IR stars (e.g. Meyer et 



al. 1998). Its formation in the circumstellar envelopes of 



OH/IR stars has been discussed in particular by Jura 
&; Morris (1985). Before the ISO mission, the far-IR ice 
bands had been observed in emission in a small number of 
sources including the OH/I R star s OH26.5, OH127.8 and 
OH231.8-h4.2 (Omont et al. |199C| and references therein). 
ISO spectra have shown the 43- and 62-mn ice bands in 
emission in various objects (Barlow 1998| ), s uch as the 
planetary nebulae CPD -56°8 032 (Cohen et al. |l999|) and 



NGC 6302 (Lim et al. |l999|) , the p ost Re d Supergiant 
source AFGL 4106 (Molster et al. |l999a| ) and Herbig 



and relative amounts, and produces the resulting spec- 
trum for optically-thin emission from the individual ma- 
terials, as well the total emission from all the materials 
(shown as the thin solid line in Fig. H) . The materials used 
to fit the OH/IR star spectra were forsterite, enstatite and 
crystalline water ice; temperatures of order 50-100 K were 
used for the fitting. The detected ice features are listed in 
Table |. 

Pyroxenes (such as enstatite; see dash-dotted line in 
Fig. ||) also show strong 43-/j,m features, therefore detec- 
tion of a 43-/im band in an observed spectrum is not suf- 
ficient evidence to demonstrate the presence of H2O ice. 
However, for temperatures ^40 K, the 43-/im peak is at 
least as prominent as the 62-/im peak, so objects which 
do not have a 43-/im feature are unlikely to contain much 
water ice (unless it is very cold). Conversely, if an object 
shows the 62-/j,m feature, H2O ice is likely to be responsi- 
ble for at least part of that object's 43-/im feature. 

We claim detections of crystalline water ice emission 
based on the presence of the bands at 43 and 62 /Ltm in 
OH127.8, OH26.5 and AFGL 5379, confirming the ten- 
tative detections for the first two sources by Omont et 



Ae/B e stars (Waters & Waelkens |1998| ; Malfait et al. |l998| , 



199£ ) , while the 43-/im band has been detected in absorp- 



tion toward the highly-embedded sources AFGL 7009 and 
IRAS 19110-H1045 (Dartois et al. |l998i Cox & Roelfsema 



1999). Both of the far-IR bands can be blended with crys- 



al. (|1990|). OH32.8 also shows a 43-/im feature (Fig. ^, 
but without observations of the 62-/im feature, we cannot 
determine if ice emission is present. 

CRL 2199 and WX Psc both seem to show broad 50- 
70 /im features, but the shape of these features does not 
resemble laboratory crystalline ice features, unlike the ob- 
served features of OH127.8, OH26.5 and GL5379. In par- 
ticular, the emissivity of ice has a minimum near 55 /xm be- 
fore reaching its peak at 62 /im (see dashed line in Fig. H). 
This structure is present in the spectra of the latter three 
sources, while the CRL 2199 and WX Psc features are 
more flat-topped, with strong emission at 55 /im and no 
real evidence of a peak at 62 jiui. The attempt to fit the 
WX Psc spectrum with ice emission (Fig. ^) illustrates 
this point. The OH 104.9 features are rather weak and ill- 
defined; there may be a weak 62-/im band, but the 43-/im 
band is replaced by a broader feature peaking near 47 /im. 

If the 50-70 /im features in CRL 2199 and WX Psc 
are not water ice, what are they? They may simply be 
instrumental artefacts: Fig. || illustrates how the spectra 
are dominated by the steep downward slope of the SED, 
and that spectral features around 60 /im are small pertur- 
bations on this overall trend. (Clearly, this statement also 
holds for the features which we believe to be real.) Ex- 
periments with adopting different LWS dark currents and 
with methods for dealing with detector memory effects 
made little difference to the spectra. The clino-pyroxene 



optical constants of Koike et al. (1993) show a broad fea- 



talline silicate emission, but examination of the shapes and 
positions of the observed bands can distinguish between 
silicate and ice emission. 

The new ISO spectra have significantly better resolu- 
tion and sensitivity than the earlier KAO data. Fig. H ture around 60 /im (see also Cohen et al. 1999), but this 
shows attempts to fit the 30-90 /im region of our 
(continuum-subtracted) spectra, using a spectral synthesis 
routine kindly provided by Dr T. Lim (personal comm.). 
This routine takes absorption (or emission) efficiencies for 
materials of interest, along with user-defined temperatures 



peaks at longer wavelengths than the crystalline water ice 
band, and so is not a likely carrier for the CRL 2199 
and WX Psc features. Similarly, the 62-/im features in 
OH127.8, OH26.5 and GL 5379 do not need an additional 
long- wavelength component, so the Koike pyroxene is not 
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Fig. 5. Fits to the spectral features in the 30-90 /im region. Noisy line: continuum-subtracted spectrum; dashed line: 
crystalline water ice; dotted line: forsterite; dash-dotted line: enstatites (sum of clino- and ortho- forms); solid line: 
sum of all components. The optical constants for ice are from Schmitt et al. (1995) and those of the silicates are from 
Jager et al. (1998). The observed 43-^m features are likely to be blends of ice and pyroxene emission. 



necessary to fit these spectra. The Jager et al. (1998) lab- 
oratory data do not reproduce the broad 60-/xm band seen 
in the Koike data, so the feature may not be real. 

In general, we find that the observed features in Fig. 
are narrower than our fits, suggesting that we have over- 
estimated the continuum level in this wavelength region, 
or that the optical constants we used do not adequately 
represent the circumstellar materials. 



The sources in which we detect the ice emission fea- 
tures are the four stars with the deepest 10-fiia silicate 
absorption, and hence presumably the highest mass-loss 
rates. These same stars also show the 3.1-fiia absorption 
band (see Fig. g). OH32.8 was too faint at 3 /^m to be 
detected by the SWS, but the S.l-fim ice absorption fea- 
ture has been detected in ground-based spectra (Roche & 
Aitken 1984). The sources with self-absorbed silicate emis- 



R.J. Sylvester et al.: 2.4 - 197 fim spectroscopy of OH/IR stars 



11 



sion features, do not appear to show ice features. OH104.9, 
which shows a relatively shallow silicate absorption, may 
show a weak 3.1-^m absorption, but it is hard to discern, 
because the spectrum is noisy at short wavelengths. 

The 6.0-//m band of water ice is signifi cantly weaker 
than the 3.1/^m band (see e.g. Moore 199£ ). The spectra 
of our most heavily-obscured sources, OH32.8, AFGL 5379 
and OH26.5, show a weak depression around 6 /ini, but 
this wavelength region is very rich in gaseous H2O lines, 
making it difficult to ascribe the observed feature to ice 
absorption. 

Ice formation requires cool temperatures and sufficient 
shielding from stellar and interstellar radiation (e.g. Whit- 
tet et al. 1988). The high densities in the (general) outflow 
that accompany large mass-loss rates may provide the re- 
quired shielding. Alternatively, enhanced densities could 
be provided by the formation of a circumstellar disk in the 
superwind phase, or by inhomogeneous mass loss, such as 
is apparent in studies of H2O and OH maser clumps (e.g. 
Richards et al. 1999). As discussed by Omont et al. ( 1990D , 
the presence of the 63-/im band requires that the water ice 
is at least partially crystalline, implying that the ice re- 
mained relatively warm (^100 K) for long enough to allow 
crystalline reorganization to take place. 

Optical depths and column densities for the de tected 
3.1-/im features are given in Table y. Meyer et al. ( 1998 ) 
have proposed that the ice column density correlates bet- 
ter with the ratio of mass-loss rate to luminosity (M/L) 
than with M alone. Adopting reasonable estimates (based 
on values in the literature) for these parameters, our re- 
sults support the relation between ice column density and 
log M/L proposed by Meyer et al (see their Fig. 3). How- 
ever, given the uncertainties in both parameters, and that 
the luminosity changes significantly with the variability 
phase, the relationship should be treated with some cau- 
tion. 

OH32.8 shows the 3.1-/im ice band, and an ll-/im ab- 
sorption feature in the wing of the silicate absorption fea- 
ture (Roche & Aitken 1984) which was attributed to the 
libration mode of water ice. Justtanont & Ticlens (1992) 
were able to model ground-based and IRAS observations 
of this source using silicate grains with water-ice mantles, 
which give a much broader 10-/im absorption feature than 
do bare silicate grains. 

The broad ll-/xm feature is clearly seen in the five 
sources with strong 10- /im absorption (Fig. m. It appears 
strongly in OH26.5, OH104.9, OH127.8 and OH32.8, and 
as an inflection near 11.5 /im in AFGL 5379. The con- 
tribution of this feature to the overall 10-/im absorption 
profile therefore does not appear to correlate fully with 
the presence of the other water ice bands: AFGL 5379 
shows strong far-IR ice emission and 3.1-/zm absorption, 
but only weak ll-/im absorption, while OH104.9 shows 
strong ll-/im absorption but has weak or absent far-IR 
and 3.1-/im features. 



Smith & Herman ( |1990| ) found an ll-/im absorp- 
tion feature in the spectrum of another OH/IR star, 
OH138.0-I-7.3, which does not show any ice absorption at 
3.1 /um. Since the 3.1-/xm stretching mode is intrinsically 
stronger than the libration mode, Smith & Herman con- 
cluded that the 11-^m feature observed towards OH138.0 
is not produced by water ice, and suggested that it is 
due to partially-crystalline silicates. Another possibility is 
that spectra like that of OH138.0 are the abso rption coun- 
terpart of the Little-Marenin & Little (|l99C| ) 'SiH-f' or 
'Broad' emission features, which show an emission com- 
ponent at ^11 /im on the wing of the silicate feature. 
These features have been ascribed to crystalline silic ates o r 
amorphous alumina grains (see e.g. Sloan & Price 1998 ). 
Clearly, full radiative-transfer modelling would be useful 
to determine whether ice mantles can indeed explain the 
range of ll-/Ltm features seen, or whether other grain com- 
ponents are necessary. 

The presence of strong water ice features in our spectra 
indicates that a substantial amount of the H2O in the cir- 
cumstellar envelopes may be depleted into the solid phase. 
This would decrease the amount of gas-phase H2O (and 
photodissociated OH) in the outer regions of the circum- 
stellar envelope which can be detected by maser and ther- 
mal emission. Water maser lines are observed to be rela- 
tively weaker in OH/IR sta rs tha n in objects with lower 
mass-loss rates (e.g. Likkel 1989 ). Collisional quenching 
due to the high densities in the inner parts of the outflow 
is thought to suppress the maser action; our results indi- 
cate that depletion into the solid (ice) phase may also play 
a role. 

CO ice shows features near 4.7 fim (e.g. Chiar et 
al. 199£): these are not seen in our spectra, but a broad 



absorption band around 4.3 /im, due to gas-phase CO is 
seen. This band is significantly broader than the 4.27- 
/im CO2 ice absorption feature seen in molecular clouds 



(e.g. de Graauw et al. 1996b). We see no evidence of CO2 
ice at 4.27 /im. CO2 ice shows another strong feature at 
15 /im; our spectra show some structure near this wave- 
length (Fig. H), but this may be an artefact of the instru- 
ment or data-reduction process. 

3.4- Silicates 

For these objects, most of the energy radiated by the cen- 
tral star is absorbed by the circumstellar dust shell and 
re-radiated as thermal emission, predominantly by amor- 
phous silicates. Amorphous silicates have strong features 
at 10 and 18 /xm. The current objects can be ordered by in- 
creasing optical depth of those features, which agrees with 
ordering by increasing wavelength of the peak of the SED. 
When the 10-/zm feature becomes optically thick, the en- 
ergy from the central star must be re-radiated at even 
longer wavelengths, therefore the peak position shifts to- 
wards 30 /im (see Fig. p. It is believed that the low mass 
loss rate Miras evolve into high mass loss rate OH/IR 
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Table 4. Crystalline silicate features in the spectra of our 
programme stars. The features are detected by close ex- 
amination of the spectra. Enstatite (pyroxene) is indicated 
by a P, forsterite (olivine) by an O. If the feature is found 
in emission, it is indicated by a Y+, and by a Y— when 
the feature is found in absorption. 



OH/IR stars presented here are the only objects known 
to exhibit crystalline silicates in absorption outside the 
8-13 /im wavelength region. For OH32.8 this was already 
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stars, while the IRAS colors indicate that the peak of the 
dust emission shifts towards longer wavelengths (van der 
Veen & Habing 1988). With increasing mass loss rate the 
characteristic density in the wind will increase, increasing 
the optical depth towards the central star. This evolution 
can be found in the oxygen-rich AGB stars in our sample. 

In Fig. g the objects are plotted in the same order as 
in Fig. y, however, the intensities are now in F^ units, 
and normalized with respect to the measured maximum 
intensity. The 10-/im feature of Mira is completely in emis- 
sion; for WX Psc and CRL 2199, the 10-/im feature is 
partially self-absorbed. For the OH/IR stars OH104.9, 
OH127.8, OH26.5, AFGL 5379 and OH32.8, the 10-/^m 
feature is completely in absorption, with optical depths 
ranging from 1.1 for OH104.9 to 3.6 for OH32.8. A de- 
tailed analysis of the amorphous silicate features will be 
presented in a future paper (Kemper et al, in preparation) . 
As the optical depth increases, some structure becomes 
apparent in the 20-45 /im region. 

In the lower mass-loss rate objects, strong emission 
features due to amorphous silicates are present, but there 
are no obvious narrow features at A > 20 /xm. For the 
redder objects, we find that the amorphous silicates fea- 
tures are (self-)absorbed, and that some structure is ap- 
parent at wavelengths A > 20 fim. These narrow features 
can be identified as crystalline silicates, both olivines and 
pyroxenes (Waters et al. 1996). The identifications are 



based on laboratory spectra of crystalline silicates (Jager 
et al. |1998| Koike et al. |1993| , Koike & Shibai |1998D and 
similar bands seen in other objects on which detailed 
studies have been performed, i.e. AFGL 4106 (Molster et 
al. |1999aD and to HD 45677 (Voors |l99g| ). The dashed 
lines in Fig. W represent the position of some important 
crystalline silicate complexes, which are listed in Table 0. 
These crystalline silicate features are found in emission at 
the longest wavelengths but sometimes in absorption at 
somewhat shorter wavelengths, for example the 23.6 /im 
olivine feature in OH32.8 and in OH26.5 (see Fig. |). The 



reported by Waters & Molster (1999) in comparison with 
AFGL 4106. The presence of the most important crys- 
talline silicate features is indicated in Table for all the 
sources in our sample. Note that those features are de- 
tected by close examination of the spectrum; not all fea- 
tures are visible in the overview figures presented in this 
paper. Detailed modelling of the wealth of crystalline sil- 
icate features shown by the individual objects is deferred 
to a future paper. 

The crystalline silicate features tend to appear in those 
sources having greater optical depth at 10 /im. However, 
the sharpness of the crystalline silicate features shows a 
large variation. The sharpness is expected to be deter- 
mined by properties of the crystalline silicates, such as 
the presence of impurities, the shape of the dust grains, 
and irregularities in the lattice structure. The spectrum 
of AFGL 5379 does not show the sharp crystalline peaks 
found in the spectra of the other OH/IR stars and in the 
spectra of AFGL 4106 and HD 45677, but the shapes of 
its features do resemble the laboratory spectra of crys- 



talline silicates (Jager et al. 1998). This suggests that 



the dust grains around these objects exhibit differences in 
the properties of the lattice structure, such as impurities, 
holes, and edge effects due to grain size. The appearance 
of the crystalline silicate emission features in the redder 
sources in our sample confirms the relationship between 
dust crystallinity and envelope colour temperature, and 
hence mass- loss rate, identified by Waters et al. ( 1996 ). 
Specifically, there seems to be a threshold value for the 
mass loss rate above which the crystalline silicate features 
appear in the spectrum. However, above this threshold 
value, the strength and width of the features seem to be 
uncorrelated to the mass loss rate. 

Fig. R presents a more detailed overview of the loca- 
tion of the crystalline silicate features. In the upper panel, 
the upper spectrum is that of OH104.9. The solid line rep- 
resents the continuum fit obtained using the method de- 
scribed above. At longer wavelengths the spectrum shows 
emission features superposed on the pseudo-continuum. 
The 18- /im absorption feature extends to A w 26 /im. How- 
ever, as Fig. clearly shows for both plots, there are crys- 
talline emission features present within this wavelength 
region, in particular, the 17.5-20 /im complex. For refer- 
ence, the continuum-subtracted spectrum of AFGL 4106 



(Molster et al. 1999a) is also plotted. The 20.6 /im emis- 
sion feature is probably related to crystalline silicates, al- 



though it is not yet identified (Molster et al. 1999a; Voors 
1999). These OH/IR star spectra are the first to show 
crystalline silicates in emission simultaneously with amor- 
phous silicates in absorption in the same wavelength re- 
gion. 

The presence of crystalline emission features in the 
spectrum of OH104.9, at wavelengths where the amor- 
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Fig. 6. Overview of the crystalline silicates. Intensities are plotted in Fi, units, and are normalized by dividing by the 
maximum intensity. The spectra are ordered by position of the peak in the determined continuum, according to Fig. nl 
and shifted for clarity. The dotted lines indicate the wavelengths where crystalline silicate features are found in two 
or more of the plotted spectra, see text. The crystalline species are indicated by a P (pyroxene) or O (olivine); see 
Table || for details. 



phous dust component is still in absorption, implies that 
the crystalline silicate dust must have a different spatial 
distribution than the amorphous silicate dust. We consider 
two possible geometries: spherical and axi-symmetrical. 

For the case of a spherically symmetric distribution, 
the crystalline dust can have a different radial distribu- 
tion and be located further out in the envelope than the 
amorphous dust. The SWS and LWS beam sizes are much 
larger than the angular size of the dust shells of these 
OH/IR stars, so the amorphous silicate absorption can 
originate from the entire dust shell, while the crystalline 
silicate emission can arise from the cool outer layers of the 
dust shell, where the material is optically thin. If the crys- 
talline silicates are located further out, we can conclude 
that the crystalline and amorphous dust has not formed 



at the same time, but that the amorphous dust annealed 
as it moved away from the star. This could imply that 
lower mass loss rate Mira variables could in principle be 
able to form crystalline material as well, but that we can- 
not detect it because the column densities are not high 
enough. However, given that higher temperatures are re- 
quired for annealing amorphous silicates into crystalline 
silicates than arc required for the formation of amorphous 
silicates themselves, it does seem unlikely that crystalline 
silicate grains could be formed further out in an outflow 
than amorphous silicate grains. One possibility is that 
some small fraction of the particles which formed in an 
outflow (perhaps the smallest particles) immediately an- 
nealed into crystalline silicates in the inner, hottest re- 
gions. The large total column density of amorphous sili- 
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Fig. 7. The upper panel shows the spectrum of OH104.9 
(upper spectrum), together with its continuum fit (solid 
line). A fit to the 18-/xm absorption feature is also plotted. 
The lower spectrum in the upper panel is the continuum 



subtracted spectrum of AFGL 4106 (Molster et al. 1999a ). 
The dashed lines indicate a correlation between crystalline 
silicate emission features in both spectra. The lower panel 
is a close up of the upper panel, for the wavelength region 
14 - 26 /im; the spectrum of AFGL 4106 is shifted upward 
for clarity. The crystalline species is indicated on the plot 
by the symbols P (pyroxene) and O (olivine). 

cate grains would lead to net absorption in the 10- and 
18-um bands and thus obscure these hot crystalline sili- 
cates. When these particles are cooled while flowing out- 
wards, the strong emission features of crystalline silicates 
at wavelengths longer than ^ 15 /xm can be seen super- 
posed on the amorphous dust features. 

An alternative scenario for the observed behaviour is 
that OH/IR stars possess a dust disk, in addition to a more 
spherically symmetric outflow. Two alternatives suggest 
themselves: 

— Crystalline silicates that have formed and moved out 
from the inner regions of the outflow produce emission 



bands that are seen superposed against the continuum 
emission and amorphous silicate features that arise 
from the disk. This would require that those OH/IR 
stars that have the largest lG-/im amorphous silicate 
optical depths are the ones whose disks are seen most 
nearly edge-on. 
— The crystalline silicates are located in the disk, while 
the amorphous silicates are located mainly in the out- 
flow and thus have a more spherically symmetric dis- 
tribution, with the amorphous silicate absorption fea- 
tures arising from optically thick lines of sight towards 
the central star. Depending on the inclination angle of 
the system, the crystalline silicate features can be seen 
in emission. When the disk is viewed face-on, the crys- 
talline silicate features would be optically thin. Radia- 
tion from most of the disk surface reaches the observer 
via lines of sight which pass through only the outer re- 
gions of the spherical dust shell, where the amorphous 
material is not optically thick. 

Two mechanisms can be invoked to explain the high 
abundance of crystalline silicates in a disk. Below the so- 
called glass temperature, only amorphous silicates con- 
dense, while crystalline silicates can form at temperatures 
greater than the glass temperature. At the high densities 
occurring in a disk, condensation of silicates will be able to 
proceed at higher temperatures than usual (Gail & Sedl- 
mayr 1999) and the temperature range in which it is pos- 
sible to condense crystalline silicates is thus broadened. In 
our current sample there is empirical evidence for a corre- 
lation between crystallinity and density. Second, if there is 
amorphous material present in the disk, this can be trans- 
formed into crystalline material by annealing. In order to 
allow the annealing process of the silicates to proceed, 
the dust-forming region should not cool too rapidly. An 
orbiting (or slowly outflowing) disk provides the required 
stability and keeps the amorphous silicates relatively close 
to the central star for a sufficiently long time for the an- 
nealing into crystalline material to occur. Both the above 
mechanisms provide circumstances that allow the forma- 
tion of crystalline silicates, which would not be the case if 
the stellar wind removes the newly-formed silicates at the 
outflow velocity, such that they rapidly cool. 

In order to study the spatial distribution of both 
dust components in more detail, spectral mapping of the 
OH/IR stars is necessary. Then we may be able to put con- 
straints on the spatial distribution, and determine the an- 
nealing time, of the crystalline dust, using travel time and 
stability arguments. Using laboratory data on annealing 
time scales, we may be able to determine physical param- 
eters such as temperature and density in the circumstellar 
dust shell, thus helping to clarify the AGB mass loss phase 
of stellar evolution. Recent high-resolution imaging of the 



dust around evolved stars, suc h as M ira (Lopez et al. [1997D 
and VY CMa (Monnier et al 



1999) 



has shown that sub- 
stantial inhomogeneities exist in the dusty outflows; our 
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spectra suggest that the dust around OH/IR stars may 
also show complex morphologies. 

For any geometry of the circumstellar shell, the 
amorphous-crystalline volume ratio could also be affected 
by grain-grain collisions. Such collisions have long been 
recognised as important processes for the evolution of 
grains in circumstellar envelopes (see e.g. Biermann & 
Harwit 198C) and in the interstellar medium (Jones et 
al. 19961 ). The shock wave driven into the two grains by the 
collision can lead to vaporization, the formation of high 
pressure phases, melting, annealing, and shatteri ng de - 
pending on the pressures involved (cf., Tielens et al. 1994 ). 
Experiments show that at a relative collision velocity of 
~1 km s~^, mechanical effects (shattering, crater forma- 
tion) become important. Thermal effects such as crystal- 
lization, involving the intergranular nucleation of new, 
strain-free grains in a previously highly deformed ma- 
trix start at somewhat higher velocities (~ 5 km s~^; 700 
kbar) and are never very pervasive. Above r^ 7 km s~^ (1 
Mbar) , melting followed by rapid quenching leads to glass 
formation (Bauer 1979; Schaal et al. 1979). Thus, these 
experiments imply that crystallization only occurs over a 
very narrow collision velocity range and is not very effi- 
cient. Of course, if the projectile/target size ratio is large, 
even high velocity impacts will lead to a small volume frac- 
tion of the target grain passing through the regime where 
recrystallization can occur when the shock wave expands. 

Given the grain velocity profile in AGB outflows 
(Habing et al. 1994), potentially crystallizing grain-grain 



collisions will be largely confined to the acceleration zone 
at the base of the flow. An important objection against an- 
nealing through grain-grain collisions is that the cooling 
time scales of dust grains are very short compared to the 
annealing time scales. Therefore the silicate d ust gra ins so- 
lidify in the amorphous state (Molster et al. 19991: ). The 
annealing process can be described as a three dimensional 
random walk diffusion process on a cubic lattice (Gail 
1998| ), which provides an estimate of the annealing time 



scale as a function of dust temperature. At a Td — 2000 K 
the annealing time scale is '^ 1 s and strongly increases 
for lower temperatures. The cooling time scale can be de- 
rived under the assumption that the power emitted by 
the dust grain is given by the Planck function, modified 
by the Planck mean of the absorption efficiency Qabs. By 
comparing the emitted power to the internal heat of the 
grain one finds for the cooling time scales TcooI ~ 10""^ s 
for Td = 2000 K and TcooI ~ 0.02 s for Td = 800 K. In 
agreement with the experiments, only a small fraction of 
crystallized material is expected due to the grain-grain 
collisional shock loading. Final assessment of the viabil- 
ity of this mechanism for the formation of crystalline sil- 
icates in AGB outflows has to await detailed modeling of 
such grain-grain collisions in circumstellar outflows (Kem- 
per et al. in preparation). Finally, we note that the spec- 
trum of /3 Pic, which is due to a dust size distribution 
that is surely coUisionally dominated, shows little evi- 



dence for crystalline silicates. Although the 10 ^m spec- 
trum of (3 Pic (Knacke et al. |1994 Aitken et al. |1993D 
may show some spectral structure resembling that of so- 
lar system comets - characteristic for silicate minerals - 
the longer wavelength bands so prominent in cometary 
spectra arc completely absent in this source (Pantin et 
al. 1999). While this may at first sight argue against this 
mechanism, the collisional velocities in this system may 
be much lower and predominantly lead to shattering. 

4. Conclusions 

We have presented the complete SWS/LWS spectra for 
seven oxygen-rich evolved stars, together with the SWS 
spectrum of an eighth source. For the OH/IR stars, which 
have optically thick dust shells, essentially all of the stel- 
lar luminosity is radiated in the wavelength region cov- 
ered by ISO. Emission features of crystalline silicates are 
seen longwards of 15 /im in the dustier objects. Some of 
these emission features lie within the 20-/xm absorption 
feature of amorphous silicate, suggesting that the crys- 
talline and amorphous components have different spatial 
distributions. The dust shells of these sources are suffi- 
ciently cool for abundant water ice to form, as indicated 
by the near-IR absorption and far-IR emission features of 
crystalline ice. 

In a future paper (Kemper et al., in preparation) we 
will present some of the further analysis and modelling re- 
quired to determine the physical conditions and processes 
which give rise to these very rich spectra. 
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